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Abstract:  We  measured  the  optical  absorptance  of  superconducting 

nanowire  single  photon  detectors.  We  found  that  200-nm-pitch,  50%-fill- 
factor  devices  had  an  average  absorptance  of  21%  for  normally-incident 
front-illumination  of  1. 55- |_im- wavelength  light  polarized  parallel  to  the 
nanowires,  and  only  10%  for  perpendicularly-polarized  light.  We  also  meas¬ 
ured  devices  with  lower  fill-factors  and  narrower  wires  that  were  five  times 
more  sensitive  to  parallel-polarized  photons  than  perpendicular-polarized 
photons.  We  developed  a  numerical  model  that  predicts  the  absorptance  of 
our  structures.  We  also  used  our  measurements,  coupled  with  measurements 
of  device  detection  efficiencies,  to  determine  the  probability  of  photon 
detection  after  an  absorption  event.  We  found  that,  remarkably,  absorbed 
parallel-polarized  photons  were  more  likely  to  result  in  detection  events 
than  perpendicular-polarized  photons,  and  we  present  a  hypothesis  that 
qualitatively  explains  this  result.  Finally,  we  also  determined  the  enhance¬ 
ment  of  device  detection  efficiency  and  absorptance  due  to  the  inclusion  of 
an  integrated  optical  cavity  over  a  range  of  wavelengths  (700-1700  nm)  on 
a  number  of  devices,  and  found  good  agreement  with  our  numerical  model. 
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1.  Introduction 

The  short  reset-time,  low  jitter,  and  broad  wavelength  response  of  superconducting  nanowire 
single  photon  detectors  (SNSPDs)  makes  them  attractive  candidates  to  replace  other  single  pho¬ 
ton  detectors,  including  avalanche  photodiodes  [1],  in  applications  such  as  free-space  optical 
communications  [2],  quantum  cryptography  networks  [3,  4],  and  quantum  computation  [5].  To 
fully  take  advantage  of  the  boost  in  speed  afforded  by  these  properties,  SNSPDs  with  high 
system  detection  efficiencies  (SDE)  are  needed.  But  SNSPDs  typically  exhibit  SDEs  of  only 
0.2-10%  for  1. 55- qm- wavelength  single-photons  [6,  7,  8].  This  limitation  is  currently  believed 
to  be  due  to  poor  coupling  efficiency  (t ]c);  however,  this  may  not  be  the  entire  story.  To  un¬ 
derstand  where  the  losses  come  from,  we  need  to  examine  the  factors  that  contribute  to  the 
SDE. 
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The  SDE  is  a  product  of  two  lumped  quantities,  the  coupling  efficiency  (17  c)  and  the  device 
detection  efficiency  (DDE).  The  quantity  r)c  encapsulates  all  loss  mechanisms  encountered 
between  the  photon  source  and  SNSPD,  and  is  defined  as  the  ratio  of  the  number  of  photons 
that  reach  the  active  area  to  the  number  of  photons  emitted  by  the  photon  source.  It  is  relatively 
straightforward  to  realize  r)c  >  0.9  through  careful  optical  design  that  minimizes  these  losses. 
On  the  other  hand,  maximizing  DDE  is  not  as  simple. 

The  DDE  is  the  probability  that  a  photon  incident  on  the  active  area  results  in  a  voltage 
pulse.  DDE  depends  on  two  quantities,  the  absorptance  A,  and  the  probability  of  resistive  state 
formation  due  to  an  absorption  event  Pr  through 


DDE  =  PRA.  (1) 

While  PR  depends  on  the  microscopic  physics  of  the  nanowires,  A  depends  only  on  the  optical 
properties  of  our  fabricated  structures  and  the  incident  field.  A  needs  to  be  maximized  to  make 
efficient  detectors;  however,  there  are  many  ways  an  incident  photon  can  remain  unabsorbed. 
For  example,  the  photon  can  pass  through  open  gaps  between  the  nanowires  or  be  transmit¬ 
ted  through  the  subwavelength  thickness  of  NbN.  In  addition,  a  low  value  for  Pr  will  further 
suppress  DDE  (and  therefore  also  SDE). 

In  this  paper,  we  experimentally  demonstrate  that  the  optical  polarization  and  changes  in  the 
geometry  (pitch  and  fill-factor)  of  SNSPDs  impact  their  absorptance  and  efficiency.  For  exam¬ 
ple,  for  devices  with  200  nm  pitch  and  50%  fill  factor,  we  found  that  21  ±  1%  of  light  incident 
from  the  front  was  absorbed  for  parallel  polarization,  while  only  10  ±  1%  was  absorbed  for  per¬ 
pendicular  polarization.  The  absorptance  was  reduced  to  14  ±  1%  and  6±  1.5%  for  parallel  and 
perpendicular  polarizations,  respectively,  for  devices  with  the  same  pitch  but  a  25%  fill-factor. 
We  describe  both  a  numerical  model  that  predicts  the  absorptance  of  our  structures,  and  ex¬ 
periments  that  directly  measure  absorptance.  We  found  that  the  experimental  data  matches  the 
predictions  of  our  model  to  within  the  uncertainties  in  our  knowledge  of  physical  and  optical 
parameters.  We  also  measured  the  DDE  for  the  same  devices  and  compared  it  to  the  measured 
absorptance.  We  discovered,  remarkably,  that  Pr  is  different  for  photons  polarized  perpendicu¬ 
lar  (_L)  to  the  nanowires  (Pr.x)  than  for  photons  polarized  parallel  to  the  nanowires  (Pr  n).  We 
propose  a  model  that  qualitatively  explains  our  observation  of  PR.|  /  Pr.x. 

Finally,  we  also  made  measurements  to  confirm  that  our  numerical  model  and  modeling 
parameters  are  accurate  over  a  range  of  wavelengths.  We  measured  the  factor  of  enhancement 
of  DDE  due  to  the  addition  of  an  integrated  optical  cavity  for  single -photon  wavelengths  of  700- 
1700  nm,  and  found  good  agreement  with  our  numerical  model.  As  a  result,  we  have  shown 
that  optical  modelling  can  be  used  to  design  future  SNSPDs  in  new  regimes  of  operation. 

This  paper  is  organized  as  follows.  In  Section  2,  we  will  describe  a  numerical  model  that 
predicts  the  absorptance  of  our  structures.  In  Section  3,  we  will  describe  the  setup  and  exper¬ 
iments  conducted  to  measure  absorptance  on  our  devices.  In  Section  4,  we  will  describe  the 
experiments  conducted  to  measure  the  wavelength  dependence  of  the  enhancement  due  to  the 
addition  of  an  optical  cavity.  We  will  not  describe  fabrication  details  or  DDE  measurement 
details  as  these  were  discussed  in  [9]  and  [10],  respectively. 

2.  Optical  model 

We  modeled  the  absorption  process  of  photons  by  an  SNSPD  as  a  plane  wave  interacting  with 
an  infinite  grating.  These  two  approximations  are  justified  because  in  our  experiments,  the 
photon  was  launched  at  normal  incidence  and  had  plane-wave  phase-fronts  when  it  reached  the 
detectors,  and  because  our  beam  spot  was  smaller  than  the  test  grating  which  will  be  described 
in  Section  3. 
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While  analytic  methods,  such  as  form-birefringence  theory,  have  been  used  in  the  past  to 
analyze  grating  structures  such  as  wire-grid  polarizers  [11,  12],  they  break  down  for  subwave¬ 
length  grating  thicknesses,  in  particular  for  electric-field  polarized  perpendicular  to  the  wires 
in  the  grating  (i.e.  in  the  direction  of  grating  periodicity) [13,  14].  Recently,  a  numerical  tech¬ 
nique  (finite-element  analysis)  was  used  to  approach  this  problem  [15];  but  because  of  either 
a  difference  in  material  parameters  or  simulation  conditions,  their  findings  were  not  consistent 
with  our  experimental  or  theoretical  results. 


(a)  device  cross-section  X 


Fig.  1 .  Infinite  wire  grid  schematized  in  (a)  is  further  reduced  via  symmetry  to  the  unit  cell 
shown  in  (b)  in  numerical  modeling  of  the  absorptance,  A.  (The  schematic  is  not  drawn  to 
scale.)  Plots  of  the  calculated  A  as  a  function  of  fill-factor,  /  and  pitch  p  are  shown  in  (c) 
and  (d)  for  parallel  (|  | )  and  perpendicular  (±)  electric  field  polarization.  Inset  in  (d)  shows 
how  the  calculated  cross-sectional  time-averaged  electric  field  magnitude,  |£cc'|  varies  with 
position  x  for  ||  and  _L  polarizations  across  a  f  =\,  p  =  200nm  structure.  The  NbN  film 
thickness  (fNbN)  was  4  nm  and  incident  electric  field  magnitude  |£o|  =  1-  The  NbN  region 
extends  from  x  =  —  50nm  to  x  =  +50nm. 

We  used  a  different  finite-element  analysis  software  (Comsol  Multiphysics  v3.2b,  EM  mod¬ 
ule)  than  the  one  used  in  [15]  to  calculate  the  absorptance  of  the  geometry  shown  in  Fig.  1(a). 
We  took  advantage  of  symmetry  by  using  the  ‘In-plane  Hybrid-Mode  Waves’  implementation 
and  defined  the  unit  cell[16,  9,  10]  shown  in  Fig.  1(b)  with  the  desired  NbN  thickness  (t  NbN), 
pitch  ( p ),  wire  width  (w),  fill-factor  (/  =  w/p),  and  a  1.55-pm-wavelength.  The  thickness 
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of  NbNjjOy  was  held  constant  at  2  nm  because  it  was  observed  by  transmission  electron  mi¬ 
croscopy  (TEM)  to  be  the  same  thickness  for  two  different  NbN  thicknesses  [17].  The  sapphire 
wafer  thickness  was  not  included  because  an  anti-reflection  coating  was  added  to  the  sapphire- 
air  interface  in  experiments.  We  applied  periodic  boundary  conditions  on  the  left  (XY)  and  right 
(X.'Y')  edges  of  the  unit  cell,  scattering  boundary  conditions  on  the  bottom  (XX ')  and  top  (YY') 
edges.  We  applied  either  an  out-of-plane  (Em),  or  in-plane  (E  \ )  electric  field  on  the  top  edge 
(front-illumination),  which  matches  the  experimental  setup  described  in  Sec.  3.  The  material 
parameters  were  defined  in  terms  of  complex  refractive  indices  collected  either  from  ellipsom- 
etry  measurements  [18]  or  the  literature.  We  used  «NbN  =  5.23  —  i 5.82,  «NbNxov  =  2.28,  and 
^sapphire  =  1-75  [19].  We  then  applied  a  Lagrange  (quadratic)  mesh  constrained  to  0.25  nm  in 
NbN  and  NbN  AOv,  and  made  to  be  denser  at  edge  XY  than  X'Y'  to  maintain  accuracy  when 
using  periodic  boundary  conditions.  Moving  to  a  mesh  with  elements  that  were  an  order  of 
magnitude  smaller  did  not  improve  accuracy.  We  used  a  direct  linear  solver  (UMFPACK)  to 
solve  for  the  spatial  distribution  of  electric  ( E )  field,  E(x.y). 

The  absorptance  A  is  proportional  to  the  integral  of  the  time-averaged  electric  field  intensity 
in  the  NbN  film.  Because  the  electric  field  does  not  vary  significantly  over  the  thickness  of 
the  NbN,  we  can  express  A  in  terms  of  the  cross-sectional  electric  field  |  Ecci  (x)  |  in  NbN  (as 
pictured  in  Fig.  1(b))  through 

rp/2  E NbN  /  rp/2 

A=  /  Q(x,y)&x&y  /  70dx  (2) 

J-p/ 2  JO  /  J~p/2 

where  Q(x,y)  is  the  time-average  resistive  dissipation  in  the  nanowire,  and  Ia  is  the  time-average 
incident  Poynting  power  density  given  by 

/o=^(£o/Mo)1/2|Fo|2  (3) 

where  E0  is  the  time-averaged  incident  electric  field  magnitude,  eD  is  the  permittivity  of  air, 
and  fi o  is  the  permeability  of  air.  For  our  wires  which  have  typical  f^N  °f  4-6  nm,  Q(x.y)  « 
Q(x )  and  is  given  by  Ohm’s  law  to  be 

Q(x)  =  ^a>  Im[e] \ECC,  (x)  | 2 .  (4) 

Figure  1(c)  shows  the  calculated  dependence  of  A  on  /  for  constant  values  of  p  and  electric 
field  polarizations  that  were  either  parallel  (||)  or  perpendicular  (_L)  to  the  periodicity  of  the 
nanowires,  while  Fig.  1(d)  shows  A  as  a  function  of  p.  Both  figures  show  that  A  is  invariant 
with  pitch  p  when  fill-factor  /  is  kept  constant,  and  that  A  can  be  maximized  by  narrowing  the 
gaps  between  nanowires.  We  understand  that  A  |  remains  constant  with  varying  p  because  the 

electric  field  intensity  is  continuous  across  the  air-NbN  boundaries  as  a  result  of  tangential  E..i:, 

continuity,  as  shown  in  Fig.  1(d)  inset.  Only  the  fraction  of  the  electric  field  magnitude  |£^,| 
in  NbN  compared  to  air  (i.e.  /)  changes  Ay.  On  the  other  hand,  decreasing  p  while  keeping 
/  constant  decreases  A  |  .  This  effect  is  observed  because  there  are  more  air-NbN  interfaces 
per  unit  area  for  a  grating  with  a  smaller  pitch  than  one  with  a  larger  pitch,  and  the  boundary 
conditions  dictate  a  lower  /:(.(,  in  the  NbN  at  each  air-NbN  interface. 

The  plots  also  illustrate  an  important  limitation  that  A  poses  to  the  photon  detection  process: 
the  maximum  A  (and  therefore  also  DDE)  can  not  exceed  30%  with  the  configuration  shown. 
There  are  two  ways  that  A  can  be  increased  that  have  already  been  demonstrated:  (1)  by  using 
back-illumination  (i.e.  through  sapphire  edge  XX')  instead  of  front-illumination  (i.e.  through 
air  edge  YY'),  and  thereby  reducing  the  index  mismatch  with  NbN,  A  can  be  increased  up  to 
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45%;  (2)  by  fabricating  an  optical  cavity  designed  to  intensify  the  field  in  the  NbN  nanowires 
[10].  An  understanding  of  the  impact  of  SNSPD  geometry  and  optical  polarization  on  A  and  P r 
will  be  helpful  in  finding  other  ways  to  improve  A. 

We  will  now  describe  how  we  determined  the  optical  absorptance  of  the  devices  and  how  the 
results  compared  to  the  model.  We  will  also  compare  our  measurements  of  the  absorptance  to 
the  DDE  measurements  of  the  same  devices. 

3.  Absorptance  measurements 

Figure  2  shows  a  schematic  of  the  optical  apparatus  we  used  to  measure  absorptance.  Free-space 
optics  were  attached  to  a  three-axis,  computer-controlled  stage  (motorized  actuator:  Newport 
FTA-HS,  motion  controller:  Newport  ESP300).  The  stage  could  be  affixed  to  a  cryostat  spe¬ 
cially  designed  to  allow  free-space  optical  coupling  to  a  cold  sample,  or  be  used  without  a 
cryostat  for  room-temperature  measurements.  We  mounted  a  single-mode  optical  fiber  (Corn¬ 
ing  SMF-28)  to  the  stage  that  was  illuminated  by  a  1.55-  pm - w a ve length  diode  laser  (Thorlabs 
S1FC1550)  and  polarized  using  a  state  of  polarization  (SOP)  locker  (Thorlabs  PF100S).  The 
SOP  locker  was  preferable  to  a  manual  polarizer  because  it  provided  mechanical  isolation  to 
our  setup  when  we  varied  polarization  and  permitted  automation  of  the  data  acquisition  process. 
The  polarized  light  was  collimated  and  split  equally  into  two  arms.  In  one  arm,  a  Ge-based  sen¬ 
sor  (Thorlabs  S122A)  connected  to  a  power  meter  (Thorlabs  PM  100)  was  used  to  measure  the 
incident  light  intensity,  mainly  for  troubleshooting  purposes.  In  the  other  arm,  a  long-working- 
distance  microscope  objective  (Mitutoyo  M  Plan  APO  NIR  20X)  was  used  to  focus  the  light 
onto  a  device.  Fight  reflected  from  the  device  was  collected  and  measured  using  another  power 
meter.  The  transmitted  light  was  measured  by  a  third  power  meter.  We  controlled  the  power 
meters,  SOP  locker,  and  motion  controller  using  the  Instrument  Control  Toolbox  (ver.  2.4)  for 
Matlab  (ver.  7.2/R2006a)  and  custom  software. 


optical  power  meter 
(transmitted  light) 


Fig.  2.  Schematic  of  the  optical  setup  used  to  measure  the  optical  absorptance.  A  state  of 
polarization  (SOP)  locker  was  used  to  set  the  input  polarization  of  the  incident  1.55-pm- 
wavelength  light.  Light  was  focused  to  a  4  ±  0.5  pm  (FWHM)  spot  onto  a  device  and  the 
incident,  reflected,  and  transmitted  power  were  measured.  The  absorptance  was  calculated 
from  these  measurements  and  earlier  calibrations  discussed  in  the  main  text. 

It  was  necessary  to  fabricate  devices  which  were  specially  designed  to  facilitate  both  ab- 
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sorptance  and  DDE  measurements.  A  typical  device  is  pictured  in  Figure  3  where  the  ac¬ 
tive  area  (3  pm  x  3  pm)  of  the  SNSPD  is  centrally  located  within  a  large  grating  structure 
(30  pm  x  10  pm).  The  small  active  area  facilitated  DDE  measurements  because  it  enabled  bet¬ 
ter  uniformity  in  electrical  response  across  the  device.  The  larger  grating  structure,  which  has 
the  same  p  and  /  as  the  smaller  active  area,  facilitated  measurements  of  A  because  it  closely 
matched  our  simulation  geometry.  Thus,  five  groups  of  devices  with  identically  sized  large  grat¬ 
ings  and  small  active  areas  were  fabricated  with  different  (/,  p  [nm])  combinations  of  (  |,200), 
( j,300),  (^, 400),  ( j,  150),  and  (^,200)  on  the  same  chip. 


Fig.  3.  Scanning  electron  micrograph  of  a  50%  fill-factor,  200-nm-pitch  SNSPD  is  shown 
in  (a).  The  central,  photon-detecting  region  of  the  SNSPD  structure  outlined  in  (a)  is  mag¬ 
nified  in  (b).  The  extended  parallel  grating  structure  outside  the  active  area  was  necessary 
for  proximity-effect  correction  in  the  fabrication  process,  but  was  also  useful  in  expanding 
the  optically  testable  region  of  the  device. 

The  procedure  for  collecting  absorptance  data  was  as  follows.  We  first  used  the  imaging 
capability  of  our  system  to  determine  the  working  distance  used  in  our  experiments.  To  find 
the  working  distance,  we  scanned  the  laser  across  a  sapphire-gold  edge  and  varied  working 
distances  between  scans.  We  measured  the  reflected  and  transmitted  light  intensity  as  a  function 
of  position  for  each  working  distance.  We  fitted  the  measured  intensity  profiles  to  theoretical 
profiles  that  were  calculated  by  convolving  a  unit  step  function  with  a  Gaussian  function,  and 
were  thus  able  to  estimate  the  full-width  half-max  (FWHM)  of  the  beam  profile  at  different 
working  distances.  The  minimum  FWHM  for  our  setup  was  4  ±0.5  pm.  By  operating  at  a 
working  distance  that  corresponded  to  the  minimum  FWHM,  we  ensured  that  98%  of  the  laser 
intensity  was  incident  on  the  large  grating  structure. 

We  then  measured  the  reflected  and  transmitted  power  from  gold  and  sapphire  in  the  vicinity 
of  the  device.  These  measurements  were  repeated  over  polarizations  spaced  by  10°  on  the 
Poincare  sphere.  We  calibrated  the  measured  reflected  and  transmitted  power  to  the  theoretical 
reflectance  and  transmittance  ( R ,  T)  of  gold  [20,  21,  22]  (0.98,  0)  and  sapphire  (0.07,  0.93) 
calculated  from  their  refractive  index  values  at  A  =  1 .55  pm.  This  calibration  procedure  corrects 
for  absorption  or  reflection  from  optical  elements  (e.g.  optical  windows).  We  then  measured  the 
reflected  and  transmitted  power  with  the  laser  centered  in  the  active  area  of  the  SNSPD,  and 
calculated  R,  T,  and  A  =  1  —  R  —  T.  We  then  repeated  our  measurements  on  10-13  identical 
devices  in  each  of  the  five  groups  of  devices. 
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Figure  4  shows  our  measurements  of  An  and  Aj_  at  room  temperature.  Room  temperature 
measurements  did  not  differ  significantly  from  measurements  of  A  made  at  6  K,  where  the  re¬ 
flectance  and  transmittance  were  calibrated  using  the  room  temperature  refractive  index  values 
of  sapphire  and  gold  at  A  =  1 .55  |um.  We  attempted  to  fit  the  model  by  using  w  and  t  NbN  as  free 
parameters,  and  assuming  ?NbN  was  between  4  nm  and  6  nm.  Our  fit  required  wire  widths  that 
were  systematically  10-15%  larger  than  the  nominal  w  used  in  our  electron-beam-lithography 
pattern.  There  are  many  uncertainties  in  our  model  parameters  that  can  contribute  to  this.  One 
cause  is  that  we  may  have  inadvertently  exposed  the  resist  adjacent  to  the  intended  struc¬ 
ture  causing  w  to  be  larger  than  desired.  We  conducted  scanning  electron  microscopy  of  our 
nanowires  and  found  that  w  predicted  by  the  model  were  within  the  uncertainty  caused  by  in¬ 
creased  secondary  electron  emission  from  the  edges  of  the  nanowires.  However,  there  are  other 
factors  that  can  contribute  to  a  fitting  error.  One  factor  is  that  the  NbN  refractive  index  we  used 
in  our  simulations  may  not  be  accurate  because  the  refractive  index  was  measured  on  a  thicker 
film  (12  nm)  than  the  film  we  used  to  fabricated  our  devices  on  (4-6  nm).  In  addition,  the  films 
were  from  different  growth  batches  so  may  have  had  slight  differences  in  optical  properties. 
Another  factor  is  that  our  simple  model  of  infinite  gratings  may  not  be  sufficient,  i.e.  a  3D 
model  may  be  needed.  Considering  these  uncertainties,  the  qualitative  and  quantitative  fit  of 
our  measurements  with  the  model  was  good. 


w=100  nm  w=100  nm  w=100  nm  w=50  nm  w=50  nm 
p= 200  nm  p= 300  nm  p=400  nm  p=150  nm  p= 200  nm 


nominal  dimensions 


Fig.  4.  Statistical  plot  of  the  measured  parallel  (||)  and  perpendicular  (±)  absorptance  for 
devices  with  different  pitch  p  and  fill-factor  /.  Each  data  symbol  represents  measurements 
of  10-13  devices  for  the  polarization  that  yielded  the  maximum  (||)  or  minimum  (±)  ab¬ 
sorptance.  The  arrows  indicate  the  calculated  absorptance  for  structures  with  fitted  wire 
widths  (from  left  to  right)  w=104  nm,  108  nm,  1 14  nm,  58  nm,  and  55  nm,  NbN  thickness 
tNbN=4.35  nm,  and  nominal  values  of  pitch  p. 
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3.1.  Comparison  of  device  detection  efficiency  to  absorptance 

A  comparison  of  device  detection  efficiency  to  absorptance  can  answer  one  important  question 
about  our  device:  is  every  absorbed  photon  actually  detected?  To  answer  this  question,  we  meas¬ 
ured  both  DDE  (using  the  measurement  apparatus  described  in  [10]  but  with  front-illumination 
instead  of  back-illumination)  and  A  for  the  same  devices.  The  two  quantities  are  plotted  against 
each  other  in  Fig.  5  where,  for  reference,  we  have  also  plotted  lines  with  constant  slope  P r.  It 
should  be  noted  that  the  value  for  DDE  is  measured  with  a  relative  accuracy  of  ±25%  due  to 
uncertainties  in  fiber  output  power  calibration.  Some  features  of  this  plot  were  as  expected;  for 
example,  all  of  the  devices  had  higher  absorptance  than  device  detection  efficiency  U\  <  1). 
But  there  were  other  features  that  were  unexpected.  We  will  now  discuss  these  other  features 
and  see  how  they  yield  an  unexpected  new  insight  into  the  microscopic  physics  of  SNSPDs. 
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Fig.  5.  Plot  of  device  detection  efficiency  as  a  function  of  absorptance  for  the  same  devices. 
Dotted  lines  have  a  constant  slope  given  by  t h,  the  probability  of  resistive  state  formation. 


There  are  two  remarkable  features  in  the  data  plotted  in  Fig.  5.  First,  If  had  a  peak  value  of 
0.85  21  t°r  parallel  polarization.  The  natural  question  one  can  ask  is  what  happens  to  absorbed 

photons  that  do  not  lead  to  detection  events?  One  possible  explanation  is  that  the  devices  were 
in  some  way  not  biased  at  their  ’’true”  critical  currents,  i.e.  that  they  were  constricted  in  some 
way  so  that  the  absorption  of  photons  does  not  cause  a  voltage  pulse.  Another  explanation  is 
that  there  is  some  other  absorbing  medium  that  was  not  accounted  for.  We  aim  to  address  this 
issue  in  future  experiments. 

A  second  remarkable  feature  of  the  data  in  Fig.  5  is  that  If  for  parallel  polarization 
(0.55  ±8;J2  <  PR  ||  <  0.85  I,];] | )  is  different  than  PR  for  perpendicular  polarization  (0.35  Iggg  < 
Pr.±  <  0.75  1  [];];] ).  While  we  do  not  yet  completely  understand  the  origin  of  the  disparity  be¬ 
tween  If  and  If  | ,  we  have  developed  a  model  that  can  form  a  starting  point  for  further 
analysis. 

The  main  assumption  of  this  model  is  a  position-dependent  (but  polarization-independent) 
probability  density  of  resistive  state  formation,  where  x  is  the  distance  across  the 

nanowire.  The  second  element  of  this  model  is  the  time-averaged  resistive  dissipation  Q(x ) 
which  we  know  to  be  both  position-dependent  and  polarization-dependent.  These  two  elements 
can  be  combined  to  give  the  device  detection  efficiency  and  absorptance  for  parallel  and  per- 
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pendicular  polarization: 


We  can  use  this  equation,  Eq.  (1),  and  a  calculated  Q{x)  to  find  a  y/(x)  that  gives  a  different 
Pr  for  parallel  and  perpendicular  polarizations.  Using  our  data,  we  found  that  y/(x)  needs  to  be 
larger  at  the  edges  of  the  nanowire  than  in  the  center  to  explain  PR  i  >  PK  \  ,  however  the  five 
geometries  we  investigated  did  not  provide  enough  resolution  to  fit  y/(x)  to  a  specific  shape  at 
this  point. 

4.  Verification  of  the  model  through  measurement  of  the  cavity  enhancement  factor 

In  the  previous  section,  we  verified  that  our  model  predicted  the  absorptance  for  our  structures 
to  within  experimental  uncertainties.  But  this  verification  was  done  only  at  a  single  wavelength 
of  1550  nm,  mainly  because  all  of  our  optical  components,  i.e.,  the  laser,  fiber,  lenses,  SOP 
locker,  microscope  objective,  beam  splitter,  collimator,  and  photodiodes  were  optimized  for 
that  wavelength.  In  order  to  verify  our  modeling  parameters  for  a  range  of  wavelengths,  we 
measured  the  cavity  enhancement  factor  £  for  SNSPDs  over  a  range  of  wavelengths.  The  cavity 
enhancement  factor  £  is  the  ratio  of  the  absorptance  of  an  SNSPD  with  an  integrated  cavity  to 
the  absorptance  without  a  cavity  and  depends  on  the  refractive  indices  of  the  materials  and  the 
geometry.  £  is  also  the  factor  by  which  the  intensity  of  the  electric  field  in  NbN,  and  therefore 
Q(x),  is  increased  due  to  the  presence  of  a  cavity.  In  view  of  Eq.  (5),  and  because  y/(x)  should 
not  depend  on  the  intensity,  we  can  determine  £  by  measuring  DDE  for  an  SNSPD  with  a 
cavity  and  the  DDE  of  an  SNSPD  without  a  cavity  and  taking  the  ratio  of  these  two  quantities 
for  each  wavelength. 

We  generated  radiation  spanning  the  range  from  600-1700  nm  using  a  supercontinuum 
source  (Toptica  photonics).  This  source  consisted  of  a  power-amplified,  modelocked  fiber  laser 
with  100  fs  pulses  at  1550  nm,  which  could  either  be  coupled  into  a  highly-nonlinear  step-index 
fiber  (generating  a  frequency  comb  from  ~600- 1 050  nm)  or  first  doubled  using  a  periodically- 
poled  Lithium-niobate  crystal  and  then  coupled  into  a  photonic  crystal  fiber  (Blazed  photonics) 
(generating  a  frequency  comb  from  ~  1150-1700  nm).  We  used  a  Pellin-Broca  prism  to  se¬ 
lect  a  wavelength  band  out  of  one  or  the  other  of  these  outputs  (with  a  FWHM  AA /A  ~  50) 
which  we  then  coupled  into  an  optical  fiber  (SMF-28,  single  mode  from  ~  1 200- 1 600  nm).  The 
wavelength  of  the  light  output  from  the  fiber  could  then  be  tuned  by  rotating  the  prism. 

We  used  this  tunable  output  to  measure  the  enhancement  £.  The  most  obvious  way  to  do  this 
would  simply  be  to  make  calibrated  DDE  measurements  at  all  wavelengths  before  and  after 
adding  the  cavity.  However,  we  chose  not  to  attempt  this,  since  we  had  no  way  to  accurately 
calibrate  optical  power  levels  at  or  near  the  single-photon  level  over  such  a  wide  wavelength 
range.  To  avoid  the  need  for  this  calibration,  we  instead  measured  only  detection  efficiency 
ratios  for  pairs  of  detectors  side-by-side.  To  do  this,  we  fabricated  a  chip  with  225  pairs  of 
3x3.3  qm  detectors,  where  the  two  in  each  pair  could  be  read  out  separately,  and  were  spaced 
by  only  100  qm.  We  measured  critical  currents  and  room-temperature  resistances  of  all  450 
detectors,  and  identified  ~  50  pairs  for  which  these  values  were  within  5%  of  each  other.  We 
then  measured  (calibrated)  DDEs  at  1550  nm  for  this  subset  of  detectors,  and  further  narrowed 
the  experimental  sample  to  the  pairs  having  DDEs  within  5%  of  each  other.  Then,  we  added 
optical  cavities  to  one  detector  of  each  pair,  and  re-measured  at  1550  nm.  Although  we  com¬ 
monly  observe  that  DDE  for  a  given  device  can  vary  between  cooldowns  (particularly  if  some 
additional  processing  has  occurred),  these  variations  are  nearly  always  strongly  correlated  with 
a  variation  of  the  critical  current  (a  reduced  DDE  correlates  with  a  reduced  critical  current). 
So,  we  selected  from  our  pairs  of  detectors  only  those  where  the  critical  currents  remained  the 
same  to  within  5%  after  the  addition  of  cavities,  and  where  the  detector  which  did  not  have  a 
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cavity  added  had  the  same  DDE  as  previously  to  within  5%.  Lastly,  we  further  reduced  our  ex¬ 
perimental  sample  by  selecting  only  the  subset  of  devices  which  were  relatively  unconstricted 
[23],  In  the  present  context,  this  corresponded  to  those  devices  with  DDE  at  1550  nm  >17% 
for  front  illumination.  Our  final  sample  consisted  of  five  detector  pairs.  For  these  ten  detectors, 
we  measured  count  rates  at  each  of  a  sequence  of  wavelengths  spanning  the  entire  accessible 
range.  By  measuring  both  detectors  in  a  pair  in  succession  without  adjusting  the  source  in  any 
way,  we  ensured  that  the  optical  power  incident  on  the  two  detectors  was  nearly  identical,  and 
therefore  that  their  count  rates  could  be  compared  directly  to  give  S .  This  precaution  was  im¬ 
portant  since  the  optical  power  obtainable  in  a  given  wavelength  band  was  not  repeatable,  nor 
was  the  coupling  into  the  SMF-28  optical  fiber,  which  was  not  single-mode  over  a  fraction  of 
our  wavelength  range.  Given  the  latter  issue,  we  restricted  ourselves  to  the  wavelength  range 
700-1700  nm  where  the  fiber  output  could  be  fairly  well  polarized  (extinction  ratio  >  20). 

Figure  6  shows  the  enhancement  measured  in  this  manner,  and  a  calculation  of  S  as  a  func¬ 
tion  of  wavelength.  We  can  see  that  the  data  has  good  experimental  agreement  with  the  cal¬ 
culation.  Inset  Fig.  6(a)  shows  the  geometry  that  was  used  in  the  calculation  while  Fig.  6(b) 
shows  a  plot  of  the  real  and  imaginary  parts  of  refractive  index  for  NbN,  NbN  vOy,  and  HSQ. 
Measurements  of  the  refractive  indices  of  NbN  and  NbN  xOy  were  made  at  room  temperature 
by  J.  A.  Woolam,  Inc.  using  spectroscopic  ellipsometry  on  a  12-nm-thick  NbN  film  deposited 
on  a  sapphire  wafer.  Measurements  of  refractive  index  of  HSQ  were  also  performed  by  J.  A. 
Woolam,  Inc.  using  spectroscopic  ellipsometry.  n  and  k  for  gold  were  taken  from  Ref.  [21]. 


Fig.  6.  Plot  of  enhancement  factor,  S  as  a  function  of  wavelength  A.  The  dotted  line  shows 
the  calculation  that  was  carried  out  on  the  unit  cell  shown  in  inset  (a)  using  values  for  the 
refractive  indices  shown  in  inset  (b). 


5.  Conclusion 

In  this  paper,  we  reported  measurements  of  the  absorptance  of  SNSPDs  and  showed  how  A 
changed  with  the  fill-factor  and  pitch  of  SNSPDs  and  optical  polarization.  We  confirmed  that 
our  numerical  model  made  accurate  predictions  for  a  range  of  geometry,  polarization,  and  wave¬ 
length.  We  found  that  200-nm-pitch,  50%  fill-factor  devices  had  an  average  absorptance  of 
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21  ±  1%  for  light  polarized  parallel  to  the  nanowires,  and  only  10  ±  1%  for  perpendicularly- 
polarized  light.  This  disparity  in  polarization-sensitive  absorptivity  was  even  more  evident  in 
lower  fill-factor  and  narrow  wire-width  devices,  where  we  measured  that  parallel-polarized  pho¬ 
tons  were  more  than  5  times  as  likely  to  be  absorbed  over  perpendicularly  polarized  photons. 
We  also  found  that  potentially,  some  absorbed  photons  do  not  result  in  detection  events,  and 
that  this  quantity  is  smaller  for  photons  with  an  orthogonal  polarization.  These  results  present 
new  challenges  that  need  to  be  understood  and  overcome  if  higher  efficiency  devices  will  be 
possible. 
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